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Abstract. We suggest a scheme to reconstruct the covariance matrix of a two- 
mode state using a single homodyne detector plus a polarizing beam splitter and a 
polarization rotator. It can be used to fully characterize bipartite Gaussian states and 
to extract relevant informations on generic states. 



1. Introduction 

Bipartite (entangled) states of two modes of the radiation field are the basic tool of 
quantum information processing with continuous variables PJI21IH!- 

Bipartite states can be produced by different schemes, mostly based on parametric 
processes in active nonlinear optical media. Generation schemes are either Hamiltonian 
two-mode processes, like parametric downconversion [4. or mixing of squeezed states 
jS], or conditional schemes based on the generation of multipartite states followed by 
conditional measurements [Hj. 

Besides mean values of the field operators, the most relevant quantity needed to 
characterize a bipartite state is its covariance matrix. For Gaussian states, a class 
that encompasses most of the states actually realized in quantum optical labs, the 
first two moments fully characterize the quantum state |H]. Once the covariance 
matrix is known then the entanglement of the state can be evaluated and, in turn, 
the performances of the state itself in serving as a support for quantum information 
protocols like teleportation or dense coding. 

Entanglement is generally corrupted by the interaction with the environment. 
Therefore, entangled states that are available for experiments are usually mixed states, 
and it becomes crucial to establish whether or not entanglement has survived the 
environmental noise. As a consequence, besides being of fundamental interest, a simple 
characterization technique for bipartite states is needed for experimentally check the 
accessible entanglement in a noisy channel [3110111111121121, as well as the corresponding 
state purity and nonclassicality [T4] ITo]. 
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In this paper we suggest a scheme to measure the first two moments of a bipartite 
state using repeated measurements of single-mode quadratures made with a single 
homodyne detector. This is an improvement compared to the scheme of Ref. jlHj . 
where two homodyne detectors have been employed. The scheme involves fourteen 
quadratures pertaining to five different field modes. It can be used to fully characterize 
bipartite Gaussian states or to extract relevant informations on a generic state. 

In the next Section we introduce the notation and describe how to obtain the mean 
values and the covariance matrix starting from the statistics of suitably chosen field 
quadratures. In Section |3] a possible experimental realization is described in details. 
Section 0] closes the paper with some concluding remarks. 



2. Bipartite Gaussian states and reconstruction of the covariance matrix 



Our scheme is aimed to reconstruct the first two moments of a bipartite states. This 
represents a relevant piece of information on any quantum state of two modes and 
provide the full characterization of the quantum state in the case of Gaussian signals. 
Gaussian states, i.e. states with a Gaussian characteristic function, are at the heart 
of quantum information processing with continuous variables. The basic reason is 
that the vacuum state of quantum electrodynamics is itself a Gaussian state. This 
observation, in combination with the fact that the quantum evolutions achievable with 
current technology are described by Hamiltonian operators at most bilinear in the 
quantum fields, accounts for the fact that the states commonly produced in laboratories 
are Gaussian. In fact, bilinear evolutions preserve the Gaussian character of the 
vacuum state [To] . Furthermore, recall that the operation of tracing out a mode 
from a multipartite Gaussian state preserves the Gaussian character too, and the 
same observation is valid when the evolution of a state in a standard noisy channel 
is considered. 

We denote the two modes under investigation by a and b. In the following we 
assume that a and b have equal frequencies and different polarizations. The Cartesian 
operators q^ and pk, k = a,b can be expressed in terms of the mode operators as follows 
1 , + . i . 

ga= V2 ' Pa= V2 

and analogously for q b and pi,. The covariance matrix of a two-mode state is a real 
symmetric positive matrix defined as follows 



,t 
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Ap a q b 
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Aq b p b 


\ Ap b q a 
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Ap b q b 





(2) 



where AX 2 = (X 2 ) - (X) 2 and AXY = ±([X,Y] + ) - (X)(Y) denote the variance of 
the observable X and the and mutual correlations between the observables X and Y 
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respectively. [X, Y] + — XY + YX denotes the anticommutator between the operators 
X and Y. Throughout the paper (X) will denote the ensemble average (X) = Tr [R X], 
R being the density matrix describing the two-mode state. The characteristic function 
of a quantum state R is defined as the expectation values x(Ai, A2) = (D(Xi) <g> D(X 2 )) 
where Xj e C, j — 1, 2 and D(X) = exp {Xa^ — A*a} is the displacement operator. The 
most general bipartite Gaussian state corresponds to a characteristic function of the 
form 

X (A)=exp{-iA T <7A-a T X} , (3) 

where A = (Re [Ai], Im[Ai], Re [A 2 ], Im[A 2 ]) T and (•••) T denotes transposition. The 
vector X = ((q a ), (p a ), (qb), (<lb)) T contains the mean value of the Cartesian mode 
operators. The characteristic function fully specify a quantum state, i.e. any 
expectation value may be obtained as a phase space integral. Since for a Gaussian 
state the first two moments specify the characteristic function, their knowledge fully 
characterize a bipartite Gaussian state. 



2.1. Covariance matrix from quadrature measurement 

For the sake of simplicity, we rewrite the covariance matrix as follows: 

a = V -M 

where the variance V and the mean M matrices may be written as 
/ {<&) \{[Pa,q a ]+) (q a qb) (q a Pb) 



V = 



l(\Pa,q a ]+) (Pa) (PaQb) 

(Qbqa) (qbq a ) (<j£) 

V (PbVa) (PbPa) l(\p b ,qb\+) (Pi 



(PaPb) 
\{[<lb,Pb\- 



and 



M = 



I (Qa) 2 (Pa)(q a ) (q a )(qb) (q a ){Pb) \ 

(Pa}(qa) {Pa} 2 {Pa}(qb} (Pa) (Pb) 

((lb}(qa} (qb}{q a ) (qb) 2 (qb}{pb) 
\ (Pb)(q a ) (Pb)(p a ) (Pb)(q b ) (Pb) 2 J 
Once defined the quadrature operator of the mode k, namely 

e i(t> + k e'^ 



V2 



we use the following conventions: 



x k = ^fc,0 5 Dk = %k,Tr/2 1 
Zj~ = ^Cfc j7r /4 , tfc = iCfc,— 7r/4 • 



(4) 



(5) 



(6) 



(7) 



(8a) 
(8b) 
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The matrix M only contains the first moments and can be reconstructed by measuring 
the four quadratures x^ and jjk, k = a,b. We have 



(qk) = (x k ) , (Pk) = (Vk) ■ 



(9) 



In order to reconstruct the variance matrix V more quadratures are needed. Let us 
introduce the modes 



a + b 



d 



V2 



ia + b 



f 



ia — b 

~7T 



(10) 



If a and b correspond to vertical and horizontal polarizations, then c and d are rotated 
polarization modes at ±7r/4, whereas e and / correspond to left- and right-handed 
circular polarizations. After tedious but straightforward calculations, we have: 



v = l 

2 



( 2<z 2 > (zl) - (tl) 
(4) ~ (t 2 a ) 2(yl) 

(xl)-(xl) (xj)-(xl) 

V (vl) - (yj) (y 2 c ) - (y 2 d ) 



(4)-(^> (y 2 }~(y 2 d ) 



2(4) 
(4) - (tl) 



(4) - (tl) 
2(vl) 



(ii) 



Furthermore, since 



V U = V A1 = \ ((yl) - (y))) = (yl) - 1 ((xl) + (yf)) 
V 23 = V m = \ {{x}) - (xl)) = \ {(xl) + (yl)) - (xl) 



(12a) 
(12b) 



the measurement of the quadratures pertaining to mode / is not essential. Overall, in 
our scheme, the reconstruction of the covariance matrix requires the measurement of 
at least fourteen quadratures, e.g. the following ones (of course measuring also the f- 
quadratures, being additional independent measurements, would improve the accuracy 
of the reconstruction) 



x 



a; Vai ^O) "a; 

x b , y b , z b , t b , 

X>ci He-, %di Vdi 



Notice that the number of parameters needed to characterize a bipartite Gaussian state 
is also equal to fourteen. 



3. Experimental implementations 

In Section [21 we have proved that it is possible to fully reconstruct the covariance matrix 
cr by measuring fourteen different quadratures of five field modes obtained as linear 
combination of the initial pair. Here we consider an implementation based on the bright 
continuous-wave beams generated by a seeded degenerate optical parametric amplifier 
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Figure 1. Scheme of a possible apparatus to measure the covariance matrix of 
the bipartite (entangled) state generated by a DOPA. The two modes, a (vertical 
polarization) and b (horizontal polarization), pass through a (removable) A/4 wave- 
plate and a rotator of polarization finally, a PBS reflects the vertically polarized 
component of its input toward a homodyne detector, which measures the Xk,tf, 
quadrature. See text for details. 



(DOPA) below threshold based on a type-II nonlinear crystal [Tj\. The two collinear 
beams (a and b) exiting the DOPA are orthogonally polarized and excited in a continuous 
variable bipartite entangled state. In the following we assume a as vertically polarized 
and b as horizontally polarized. 

Since the mode / is not necessary to reconstruct the covariance matrix, we do not 
consider its selection, focusing our attention on modes a, b, c, d, and e. The mode under 
scrutiny is selected by inserting suitable components on the optical path of fields a and b, 
before the homodyne detector. Modes a, b, c, and d are obtained by means of a rotator 
of polarization R$ (namely a A/2 wave-plate) and a polarizing beam splitter (PBS), 
which reflects toward the detector the vertically polarized component of the impinging 
beam. The action of the rotator R$ on the basis {|V), \H)} is given by 

R#\V) = costf|V) -sintf|#), (13a) 
R#\H) = sm$\V) + cos &\H) . (13b) 

In order to select mode e a A/4 wave- plate should be inserted just before the rotator R# 
(see Fig. QJ). The A/4 wave-plate produces a 7r/2 shift between horizontal and vertical 
polarization components, thus turning the polarization from linear into circular. 

Table [T] summarizes the settings needed to select the five modes. Overall, the 
vertically polarized mode k arriving at the detector can be expressed in terms of the 
initial modes as follows 

k = exp{i(p} cos $ a + sin $6, (14) 

where ip = tt/2 when the A/4 wave-plate is inserted, (p = otherwise. 

Once the mode k has been selected, a homodyne detector is used to measure the 
generic quadrature x^. Homodyne relies on the controlled interference between the 
quantum beam (signal) to be analyzed and a strong "classical" local oscillator (LO) 
beam of phase 0. Indeed, to access one have to suitably tune the phase 0. The 
optimization of the efficiency is provided by matching the LO mode to the mode k. The 
mode matching requires precise control of the LO frequency, spatial and polarization 
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Mode 


A/4 

7 




a 


no 





b 


no 


+7T/2 


c 


no 


+7T/4 


d 


no 


-tt/4 


e 


yes 


+7T/4 



Table 1. Setting to select the different modes k. The table refers to the elements 
depicted in Fig. The mode a is assumed to be vertically polarized and the mode b 
horizontally polarized. 

properties. Remarkably, the detected mode is always vertically polarized, thus avoiding 
any need of tuning the LO polarization. 

4. Conclusions 

A simple scheme has been suggested to reconstruct the covariance matrix of two-mode 
states of light using a single homodyne detector plus a polarizing beam splitter and 
a polarization rotator. Our scheme requires the local measurements of 14 different 
quadratures pertaining to five field mode. It can be used to fully characterize bipartite 
Gaussian states and to extract relevant informations on generic states. Finally, we 
notice that an efficient source of polarization squeezing has been recently realized |18j . 
which might be considered as a preliminary stage for the experimental realization of the 
present characterization scheme. 
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